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We report first principles calculations of the electronic structure and electron-phonon coupling in 
the non-centrosymmetric superconductor LaNiC2. These show that the material is a conventional 
electron-phonon superconductor with intermediate coupling. There are large contributions to the 
coupling by two low frequency C non-bond-stretching modes, one of which has strong Kohn anoma- 
lies. Since LaNiC2 lacks inversion symmetry, the pairing is of dominant s-wave type with some 
mixture of p-wave character. This will give exponential decay in the specific heat. 

PACS numbers: 74.70.Dd,74.25. Jb,74.25.Kc 



LaNiC2 is a member of a family of ternary nickel car- 
bides i?NiC2 (where i? is a rare earth or Y) first re- 
ported by Bodak and Marusinji These compounds form 
in the base centered orthorhombic CeNiC2-type struc- 
ture (space group Amm2) that lacks inversion symme- 
try. Studies of magnetic properties show that the i?NiC2 
phases exhibit a variety of magnetic ordering for various 
Although, Ni is often associated with magnetism, 
e.g. in oxides and in intermetallics, the magnetism in 
these compounds is due to 4/ electrons of R with almost 
no magnetic contribution from Ni. Therefore, it is per- 
haps not surprising that LaNiC2 does not show magnetic 
ordering. Interestingly, it is instead superconducting at 
Tc = 2.7 KM 

Lee et al. found non-exponential decay of specific heat 
below Tc and based on this argued that LaNiC2 is an 
unconventional superconductori^ Since pure Ni is a fer- 
romagnet and Ni compounds are often magnetic, one 
might suppose LaNiC2 to be near magnetism and there- 
fore that a pairing mechanism involving spin fluctuations 
is operative. Spin fluctuations are not pairing for or- 
der parameters that do not change sign on the Fermi 
surface (i.e. conventional s-wave). However, Pecharsky 
et al. observed the usual exponential decay of specific 
heat below consistent with the conventional BCS 
superconductivity.^ One avenue to address the inconsis- 
tency of these two experiments is by determining the 
symmetry that is broken in the superconducting sys- 
tem. Gauge symmetry is broken in a conventional BCS 
superconductor^ while other symmetries are broken in 
unconventional onesi^ Related to this, Hillier et al. re- 
cently reported their results of muon spin relaxation mea- 
surements on LaNiC2 that indicate that time reversal 
symmetry is broken in the superconducting stated Such 
time reversal symmetry breaking requires a supercon- 
ducting state with triplet character, and does not occur in 
singlet superconductors. However, this analysis is com- 
plicated by the non-centrosymmetric lattice structure of 
LaNiCa. 



There are a number of other examples of non- 
centrosymmetric superconductors including for example 
CePtSia, CelrSig, CeRhaSi, Ulr, Li2Pd:rB, RcgW and 
Y2C3.^>iaiiiJ^ii^iiiii^ These include both non-phonon 
mediated heavy Fermion metals and electron phonon su- 
perconductors, which behave very much like conventional 
centrosymmctric superconductors. Since parity is not a 
good quantum number in non-centrosymmetric systems, 
it necessitates the modification of scheme for classifica- 
tion of Cooper pairs. For example, Yanase and Sigrist— 
argue that Cooper pairs have a state that has dominant p- 
wave with some mixture of s-wave symmetry in CePtaSi. 
Considering the Tc of LaNiC2, either case is possible; it 
could be a conventional electroii-phonon superconductor 
in analogy with Y2C3 (refs. [Hill) and RegW (ref. ffi), 
in which case it would have a Cooper pair state that has 
dominant s-wave but some mixture of p-wave symmetry 
or alternatively it could be an unconventional supercon- 
ductor in analogy with the non-centrosymmetric heavy 
Fermions such as CePtaSi (refs. IqHTqI). or near ferromag- 
netic Ulr, in which case it would have dominant p or d 
wave character, which would allow a simpler explanation 
of non-cxponcntial behavior. 

Here, we report the results of first principles calcu- 
lations of electronic structure, phonon dispersion and 
electron-phonon coupling. We find that the main con- 
tributions to the electronic structure from Ni 3d orbitals 
are away from the Fermi energy and consequently that 
LaNiC2 is non-magnetic and rather far from magnetic 
instabilities that might lead to strong spin fluctuations 
thereby leading to an unconventional superconducting 
state. On the other hand, we obtain a value of electron- 
phonon coupling constant A ~ 0.52 with logarithmically 
averaged frequency ujin ^ 207 cm^^. Using simplifled 
Allen-Dynes formula, we obtain Tc ~ 3.0 K that suggests 
LaNiC2 is an intermediately coupled electron-phonon su- 
perconductor with dominant pairing that has s-wave but 
also some mixture of p-wave symmetry. 

LaNiC2 forms in a base centered orthorhombic struc- 
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ture with La on site (0.5, m,1 — u), Ni on site (0,0,0) 
and C on site {0,v,w). In this compound, the La atoms 
form trigonal prisms, which are alternately filled by Ni 
and C dimers, hence breaking the inversion symmetry. 
Another characteristic of this structm-e is the short bond 
length of C dimers, which indicates the existence of very 
stiff bond stretching modes. If these modes were respon- 
sible for superconductivity, it would result in a moder- 
ately high electron-phonon coupling and a large logarith- 
mically averaged phonon frequency. This would mean a 
very high prefactor in the McMillan equation for Tj,. Pre- 
vious tight-binding^° and density functiona l^^i^^ studies 
of other metals with C dimers show that the electronic 
structure near the Fermi level could have substantial con- 
tributions from states with C-C antibonding character. 
If this were the case for LaNiC2, there may be a strong 
dependence of Tc on the occupation of the antibonding 
states (which could be controlled by doping) and perhaps 
also a high in optimally doped samples. 

Our electronic structure calculations were performed 
within the local density approximation (LDA) using 
the general potential linearized augmented planewave 
(LAPW) method as implemented in our in-house codei^^ 
We used LAPW spheres of radius 2.0ao for La and Ni 
and 1.25ao for C. We repeated some calculations with 
the WIEN2k code^S, as a test and found no significant 
differences. We used the experimentally reported lattice 
parameters (a = 3.959 A,b = 4.564 A, c = 6.204 A)jl but 
relaxed the internal coordinates. We obtain for internal 
parameters, u = 0.3886, v — 0.5411, and w = 0.1606. 
This yields the calculated C-C distance of 1.36 A, which 
is in reasonable accord with the reported experimental 
value^ of 1.41 A. 

The phonon dispersions and electron-phonon coupling 
were calculated using linear response as implemented 
in QUANTUM ESPRESSO code.^ The linear response 
calculations were also done using experimental lattice 
parameters, using ultrasoft pseudopotentials within the 
generalized gradient approximation (GGA) of Perdew et 
al^ An 8x8x8 grid was used for the zone integration 
in the phonon calculations, while a more dense 32 x 32 x 
32 grid was used for the zone integration in the electron- 
phonon coupling calculations. The basis set cutoff for the 
wave functions was 40 Ry, while a 400 Ry cutoff was used 
for the charge density. The internal coordinates were 
again relaxed and we obtained values that agreed well 
with the values obtained from LDA calculations. The 
GGA and LDA electronic structures were very similar. 

The calculated band structure and Fermi surface of 
LaNiC2 is shown in Figs. [T] and [U respectively. These 
agree in large scale features with recent calculations done 
by Laverock et al.j^- but the details of the Fermi surface 
and band dispersions differ significantly along some di- 
rections, yielding differences in the Fermi surface. These 
presumably reflect the use of non-full potential meth- 
ods in Ref. [2^. The corresponding density of states is 
shown in Fig. [31 The band structure shows a C 2s de- 
rived band between -15.5 and -13.5 eV and another C 2s 
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FIG. 1: Calculated LDA band structure of LaNiC2 plotted 
along the path (0,0,0) ^ (0,|,0) ^ (0,i,i) ^ (0,0, §) ^ 
(0,0,0) -> (i,0,0). 




FIG. 2: (Color online) Calculated LDA Fermi surface of 
LaNiC2. The shading is by velocity. The arrows indicate 
nesting vectors that give rise to Kohn anomalies. 



derived band between -8.5 and -6.5 eV, relative to the 
Fermi energy Ep. This is followed by a manifold of three 
bands with C 2p character, associated with C dimers, 
that extend from -7.5 to -3.8 eV. This is followed by a 
manifold of 5 bands of Ni 3c? character between -3.5 and 
-0.9 eV. These account for the majority of the Ni d char- 
acter in the density of states although there is some Ni 
d contribution near Ep, similar to the borocarbides (e.g. 
LuNi2B2C), which electron phonon superconductors 
There are two bands that cross the Fermi level and they 
have mixed Ni 3(i, La bd and C 2p antibonding charac- 
ters. This electronic structure is consistent with strongly 
bonded C dimers embedded in a metallic solid with states 
near Ep having La, Ni and C antibonding character. 

The density of states at Ep is N{Ep) = 2.6 e\~^ 
on a per formula unit both spin basis. The calcu- 
lated Fermi velocities are < v^,^ >^/^=1.79xl0'^ ni/s, 
< >i/2^1.52xl05 m/s, and < >i/2=1.48xl0^ 
m/s, indicating only modest resistivity anisotropy unless 
there is a large scattering anisotropy. The density of 
states has a prominent peak with a maximum very close 
to Ep. Therefore the stoichiometric compound is ex- 
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FIG. 3: Calculated LDA density of states of LaNiC2 on a per 
formula unit both spins basis. The projection is onto the Ni 
LAPW sphere. 




0.2 0.4 0.6 O.i 

M(Hb) 



1 1.2 1.4 





1400 r- 




1200 - 




1000 - 






s 


800 - 










3 


600 - 




400 - 




200 - 




L 




FIG. 5: Calculated phonon dispersion curves of LaNiC2 plot- 
ted along the path (0,0,0) (0,i,0) (0,i,i) (0,0, i) ^ 
(0,0,0) ^ (i,0,0). 
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FIG. 6: Calculated phonon density of states and electron- 
phonon spectral function a'^F{uj) for LaNiC2. 



FIG. 4: Fixed spin moment magnetic energies for LaNiC2 



pected to have the highest Tc similar to Y2C3, although 
in the present case the peak is broader. Doping away 
from stoichiometry will move the Ep away from the peak 
and is therefore expected to lower the T^. The moder- 
ately high value of N{Ep) might also indicate magnetism 
if the states at the Fermi level are of predominantly Ni 
d character. However, as mentioned, most of the Ni d 
states are in a manifold of 5 bands that lie between - 
3.5 and -0.9 eV relative to the Fermi level. Hence, the 
high N{Ef) of LaNiC2 does not place it near magnetism. 
This is supported by our fixed spin moment calculations. 
Fig. 0] shows energies calculated for different values of 
magnetic moment relative to the nonmagnetic case. As 
may be seen the energy scale is relatively high, which 
means that the material is not near ferromagnetism. 

The calculated phonon dispersions of LaNiC2 are 
shown in Fig. O The corresponding phonon density 
of states and Eliashberg spectral function a^Fluj) are 
shown in Fig. [S] The phonon dispersions show a set of 
11 bands extending up to ~ 530 cm~^, separated by a gap 
of ^ 775 cm^^ from 1 high frequency band that lies be- 
tween 1305 and 1355 cm~^. The bands below 180 cm~^ 
show contributions from La, Ni and C, while the bands 



between 180 and 255 cm~^ are derived mainly from mo- 
tions of Ni and C atoms. Above 255 cm~^ the bands are 
dominated by motions of C atoms. The highest frequency 
band is an almost pure C-C bond stretching mode, con- 
firming the picture of strongly bonded C dimcr inferred 
from band structure. 

The eighth branch of the phonon dispersion, which 
has mainly C non-bond-stretching character, shows Kohn 
anomalies at fc - (0,0.5,0.25), k - (0,0,0.25) and 
k ~ (0.25,0,0). These can be seen in Fig. [5] along the 
Y-T, Z-F and T-X lines. The origin of the Kohn anoma- 
lies can be seen in the Fermi surface, which shows nest- 
ing with wavevectors k ~ (0,0,0.25) and k ~ (0.25,0,0) 
as indicated by two arrows in Fig. [51 A comparison of 
phonon density of states with the with the Eliashberg 
spectral function shows that the latter is enhanced rela- 
tive to the former between 220 and 400 cm^^, which is 
the region where the seventh (which also has mainly C 
non-bond-stretching character) and eighth branches are 
placed. We obtain a value of the electron-phonon cou- 
pling Aep = 0.52. As there are 12 phonon branches, the 
average electron-phonon coupling per branch is Xavg = 
0.04. The contribution due to seventh and eighth bands 
is large with values of A7 = 0.13 and As = 0.09, re- 
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spectively. We obtain for the logarithmically averaged 
frequency — 207 cm^^. Inserting these numbers into 
the simplified Allen-Dynes formula, 

, ^ fi^m r 1.04(1 + Xep) \ 

with fi* = 0.12, we obtain Tc ~ 3 K, which is in accord 
with the experimental value of Tc = 2.7 K.**'^ While this 
close agreement may be partly fortuitous (note the pa- 
rameter /i* in the Allen-Dynes formula; 0.12 is a typical 
value, but values in the range 0.1 - 0.15 are reasonable), 
it does indicate that the superconductivity can be readily 
explained in an electron-phonon framework. 

To summarize, we report full potential calculations of 
the electronic structure, phonon dispersions of LaNiC2 
and electron-phonon coupling in LaNiC2. We find that 
conventional electron-phonon coupling picture readily de- 
scribes the superconductivity of LaNiC2. With inversion 
center this would yield a pure s-wave state. However, 
because of the non-centrosymmetric structure of LaNiC2 



a small p admixture will be present. The Ni d states 
are located mainly away from the Fermi energy, and ac- 
cordingly we do not find LaNiC2 to be near magnetism. 
Therefore, we expect the contribution of spin fluctua- 
tions to the coupling, which is pair breaking for s-wave, 
to be small. The phonon-mediated Cooper pairs are in- 
termediately coupled with large contributions made by 
two low frequency C non-bond-stretching modes, one of 
which has Kohn anomalies. Owing to the lack of inver- 
sion symmetry, the pairing will have dominant s-wave 
character with some mixture of p-wave. This is expected 
to give exponential decay of specific heat below Tc- In 
particular, there will be no line nodes in the gap function. 
Also, while the Cooper pairs also have some p-wave char- 
acter, a pure electron phonon mechanism will not yield 
breaking of time reversal symmetry due to the removal 
of centrosymmetry. 
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